arXiv: 1506.02841 vl [cond-mat.mtrl-sci] 9Jun2015 


Computational 2D Materials Database: Electronic 
Structure of Transition-Metal Dichalcogenides and 

Oxides 


This document is the unedited Author’s version of a Submitted Work that was subsequently 
accepted for publication in The Journal of Physical Chemistry C, copyright (c) American Chemical 
Society after peer review. To access the final edited and published work see 
http://dx.doi.org/10.1021/acs.jpcc.5b02950 

Filip A. Rasmussen and Kristian S. Thygesen* 

Center for Atomic-scale Materials Design (CAMD) and Center for Nanostructured Graphene 
(CNG), Department of Physics, Technical University of Denmark 

E-mail: thygesen(a)fysik.dtu.dk 


Abstract 

We present a comprehensive first-principles study 
of the electronic structure of 51 semiconducting 
monolayer transition metal dichalcogenides and - 
oxides in the 2H and IT hexagonal phases. The 
quasiparticle (QP) band structures with spin-orbit 
coupling are calculated in the GoWq approxima¬ 
tion and comparison is made with different den¬ 
sity functional theory (DFT) descriptions. Pitfalls 
related to the convergence of GW calculations for 
2D materials are discussed together with possi¬ 
ble solutions. The monolayer band edge posi¬ 
tions relative to vacuum are used to estimate the 
band alignment at various heterostructure inter¬ 
faces. The sensitivity of the band structures to 
the in-plane lattice constant is analysed and ratio¬ 
nalized in terms of the electronic structure. Fi¬ 
nally, the ^-dependent dielectric functions and ef¬ 
fective electron/hole masses are obtained from the 
QP band structure and used as input to a 2D hydro- 
genic model to estimate exciton binding energies. 
Throughout the paper we focus on trends and cor¬ 
relations in the electronic structure rather than de¬ 
tailed analysis of specific materials. All the com¬ 
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puted data is available in an open database. 

1 Introduction 

Atomically thin two-dimensional (2D) materials, 
such as graphene, hexagonal boron-nitride, and 
the transition metal dichalcogenides (TMDs) are 
presently being intensively researched due to their 
unique opto-electronic properties. The TMDs with 
the chemical formula MX 2 (X=S, Se, Te and 
M=transition metal) represent a particularly in¬ 
teresting class of 2D materials comprising both 
semi-conductors and metals.^For example, M 0 S 2 , 
MoSe 2 , WS 2 , and WSe 2 were shown to undergo 
a transition from indirect to direct band gap ma¬ 
terials when their thickness is thinned down to a 
single layerTogether with their strong interac¬ 
tion with light l^ a nd relatively high charge car¬ 
rier mobilities,i21i21 this has opened up for the use 
of few-layer TMDs in a range of applications in¬ 
cluding ultrathin field effect transistors,^ photo 
detec tors light-emitting diodes,!^ and solar 
cells .[nUU Furthermore, the lack of inversion sym¬ 
metry in the 2H monolayer structures leads to 
a spin-orbit driven splitting of the valence band 
which in turn allows for valley selective excita- 
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tion of charge carriers Adding to this the 

possibility of tuning the electronic properties by 
stra in,!^ dielectric screening ,^231 electrostatic gat¬ 
ing, nanostructuring, ^ or by combining in¬ 
dividual 2D mate rials into van der Waals het¬ 
erostructures it is clear that monolayer TMDs 
hold great potential both as a platform for funda¬ 
mental physics and as building blocks for nano¬ 
scale device applications. 

Until now, opto-electronic research in mono- 
layer TMDs has mainly focused on the Mo and 
W based compounds which have (o ptical) band 
gaps in the range 1.6-2.0 significantly 

larger than the ideal values for both photovoltaics 
and transistor applications.^ Furthermore, in or¬ 
der to advance the usage of 2D materials from the 
level of fundamental research to real applications, 
it is essential to enlarge the space of available “2D 
building blocks” beyond the handful of presently 
considered materials. To this end, not only the 
band gaps but also the absolute band edge posi¬ 
tions, effective masses, and dielectric function will 
be of key importance for predicting the usefulness 
of a given 2D material. 

The fact that the interlayer bonding in bulk 
TMDs is of very similar strength (around 
20meV/A^)l^ indicate that exfoliation of sin¬ 
gle layers should be feasible for many different 
TMDs. Indeed, liquid exfoliation of nanosheets 
of TaSe 2 , NbSe 2 , NiTe 2 , and MoTe 2 has already 
been demonstrated.^ In this regard, it is interest¬ 
ing to note that more than 40 TMDs are already 
known in the bulk form and could form the basis 
for new 2D materials.^ The stability of such 2D 
monolayers under ambient conditions is a critical 
issue, but it could be alleviated by encapsulation 
in protecting layers as recently demonstrated for 
M 0 S 2 in hexagonal boron-nitride.^^ 

In a previous work, Ataca et al. performed 
an extensive stability analysis of 88 monolayer 
TMDs and TMOs using density functional theory 
(DFT) in the local density approximation (LDA), 
and identified 52 stable compounds including both 
metals and semiconductors.^ While stability was 
their main focus, they also calculated the LDA 
band structures of the stable compounds and a 
few selected compounds using the GWq approx¬ 
imation. They concluded surprisingly, that the 
LDA provides good agreement with existing ex¬ 


periments while GWq significantly overestimates 
the band gap. This false conclusion is based on 
the common confusion between the optical and 
the quasiparticle (QP) band gaps. The former is 
probed in optical experiments and is lower than 
the QP gap by the exciton binding energy. It is 
one of the characteristic features of the atomically 
thin semiconductors that exciton binding energies 
are very large (on the order of 1 eV). This leads to 
pronounced differences between the QP and op¬ 
tical spectra both of which are well reproduced 
by many-body calculations applying the GW ap¬ 
proximation and Bethe-Salpeter equation, respec- 
tively.^^^®^ 



Figure 1: (a) Crystal structure of the monolayer 
transition metal oxides and dichalcogenides in the 
2H and IT phases, respectively. Here a denotes the 
in-plane hexagonal unit cell lattice constant and 
h is the vertical distance between the oxygen or 
chalcogen atoms, (b) In-plane brillouin zone of 
the hexagonal unit cell with high symmetry points 
and other k-points indicated, (c) Periodic table 
of the elements with the metals considered in this 
study highlighted in blue and oxygen and chalco- 
gens highlighted in yellow. 

In this paper we present an extensive first- 
principles study of the electronic structure of a va¬ 
riety of monolayer TMDs and TMOs in the 2H and 
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IT structures based on 27 different metals. For 
referenee the atomie struetures of the 2H and IT 
phases are shown in Fig. the eorresponding 
Brillouin zone with the speeial ^-points is shown 
in Fig. [T^ and the elements eonsidered are high¬ 
lighted in the periodie table shown in Fig. Out 
of 216 investigated eompounds we find 171 to be 
stable (defined by a negative heat of formation rel¬ 
ative to the standard states). These results repre¬ 
sent a eonsistent extension of the LDA-based sta¬ 
bility analysis of Ref. |33 Out of the 171 sta¬ 
ble monolayers, we here foeus on the 51 eom¬ 
pounds that are found to be non-magnetie and non- 
metallie. For these materials, we ealeulate the 
band struetures as well as the absolute position of 
the valenee and eonduetion bands edges relative to 
vaeuum using the GqWo approximation with spin- 
orbit eoupling ineluded. Convergenee of the ab¬ 
solute GoWo quasi-partiele energies is found to be 
partieularly demanding and we therefore diseuss 
this issue in some detail. The GqWo band gaps and 
band edge positions are eompared to Kohn-Sham 
DFT using different exehange-eorrelation fune- 
tionals. We find that the band gap is generally well 
reprodueed by the GLLB-SC funetional while the 
LDA provides a surprisingly good deseription of 
the band gap eenter. In eontrast, an empirieal for¬ 
mula for estimating the band edge positions from 
the electro-negativities of the eonstituent atoms, 
is found to deviate signifieantly from the first- 
prineiples results due to eharge transfer from the 
metal to the oxygen/ehaleogen atoms and asso- 
eiated eleetrostatie potential that lowers the elee- 
tronie band energies relative to vaeuum. We fur¬ 
thermore ealeulate the (statie) (^r-dependent dielee- 
trie funetion of all the eompounds and diseuss 
some basie properties of dieleetrie sereening in 
quasi 2D. The effeetive eharge earner masses are 
derived from the GqWo band struetures and used, 
together with the dieleetrie funetions, as input to 
an effeetive 2D model for the exeiton binding en¬ 
ergies. 

Overall, our results reveal a large degree of vari¬ 
ation in the eleetronie properties of the investi¬ 
gated materials. For example, the materials MX 2 
(X=S, Se, Te and M=Cr, Mo, W) have direet QP 
band gaps in the range 0.9-2.5 eV while all other 
eompounds have indireet gaps in the range 0.5- 
7.0 eV. The band gap eenters (relative to vaeuum) 


span from — 8eV for some of the oxides to above 
—5eV for the selenides and tellurides. The effee¬ 
tive masses vary by almost two orders of magni¬ 
tude as do the (^-dependent dieleetrie funetions. 

All the eomputed data ineluding relaxed strue¬ 
tures, DFT and GqWo band struetures, absolute 
band edge positions, effeetive masses, and dielee¬ 
trie funetions, are available online in the Com¬ 
putational Materials Repository (http : //cmr . 
f ysik . dtu . dk/). 


2 Computational methods 

All ealeulations were performed using the projee- 
tor augmented wave method as implemented in 
the GPAW eode.® GPAW supports three different 
types of basis sets, namely real spaee grids, nu- 
merieal atomie orbitals, and plane waves. We have 
used the latter in the present work sinee exeited 
state ealeulations with GPAW are implemented 
only for plane waves. 

2.1 Atomic structure 

The lattiee eonstants (see Figure Bi) of the 
216 monolayer TMDs and TMOs were deter¬ 
mined by a strueture relaxation using the Perdew- 
Burke-Emzerhof (PBE)I^ exehange-eorrelation 
(xe) funetional with 750 eV plane wave eut-off, 
18 X 18 X 1 Monkhorst-Pack k-point sampling, 
and 20A between periodieally repeated layers. 
Eor both the 2H and IT phases the lattiee eonstant 
of the minimal unit cell and the vertieal positions 
of the oxygen/ehaleogen atoms where relaxed un¬ 
til all forees were below 0.01 eV/A. We used the 
minimal unit eell and did not investigate sym¬ 
metry redueing distortions whieh have reeently 
been found to oeeur for some IT metallie eom¬ 
pounds.®® Sinee these distortions are driven by 
a metal to insulator transition (Peierls distortion) 
we do not expeet them to be important for the 
semi-eondueting materials whieh are the foeus of 
the present work. 

2.2 Electronic structure 

The Kohn-Sham band struetures of all eompounds 
was ealeulated self-eonsistently using the EDA, 
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PBE and GLLB-SCP xc-functionals. Spin-orbit 
coupling was included in a non self-consistent 
manner by diagonalising the total Hamiltonian 
consisting of the spin-orbit interaction (which is 
applied inside the PAW spheres) and the self- 
consistently determined Kohn-Sham Hamiltonian. 
We have found that the spin-orbit corrections 
to the band structure were unchanged (less than 
0.02 eV) if we use GqWo energies instead of LDA 
energies in the Kohn-Sham Hamiltonian. For 
ten representative materials we benchmarked the 
spin-orbit corrected LDA band structures obtained 
with GPAW against the all-electron linearised 
augmented-plane wave ELK code^^ and found ex¬ 
cellent agreement (difference within 0.02 eV) 

The QP band structures were calculated in the 
Go Wo approximation as implemented in GPAW.I^ 
We used LDA wave functions obtained from an 
exact diagonalization of the Kohn-Sham Hamil¬ 
tonian with a plane wave cut-off of 600 eV and 
30 X 30 k-points, as input for the GqWo calcula¬ 
tions. The plane wave cut-off used to construct 
the screened interaction and self-energy was var¬ 
ied between 150eV and 500 eV and extrapolated 
to infinite cut-off energy as described below. For 
all calculations, the number of unoccupied orbitals 
used to construct the screened interaction and GW 
self-energy was set equal to the number of plane 
waves. The frequency dependence was repre¬ 
sented on a non-linear grid from OeV to the en¬ 
ergy of the highest transition included in the ba¬ 
sis with a gradually increasing grid spacing start¬ 
ing at 0.1 eV and reaching 0.2eV at CO = 15eV. 
The frequency grid typically contained 300 to 350 
grid points. The PAW potentials applied in this 
work include semi-core states, i.e. atomic states 
down to at least 1 Hartree below vacuum, while 
deeper lying states are included in the frozen core. 
The frozen core states are included in the exchange 
contribution to the GW corrections. 

In Ref. |43] we demonstrated the importance of 
using a truncated Coulomb interaction in GW cal¬ 
culations of 2D materials. In the present study we 
have used the Wigner-Seitz truncation scheme.^ 
For a representative set of materials we have 
checked that the QP band gaps (and more gen¬ 
erally the absolute band edge positions) changes 
by less than 0.1 eV when the k-point grid is in¬ 
creased from 30 x 30 to 45 x 45. We note in pass¬ 


ing that most previous GW calculations for 2D 
syste ms have applied significantly smaller k-point 
grids. 1361321521 explained in Ref. |43l the physi¬ 
cal reason for the slow convergence with k-points 
is the strong (^r-dependence of the dielectric func¬ 
tion of a 2D semiconductor: While £{q) for a 3D 
semiconductor tends smoothly to constant value 
for ^ 0, e(^) = 1 -1- 0{q) for a 2D system, see 

Figure As a consequence a denser k-point grid 
is required to capture the variation in £{q) around 
q = 0. For example, the GqWo band gap of 2H- 
M 0 S 2 is reduced by 0.4 eV when increasing the 
k-point grid from 15 x 15 to 30 x 30. Since the 
strong variation in £{q) is limited to a small re¬ 
gion around (? = 0, it is sufficient to sample the 
screened interaction W(q) on a fine grid in this re¬ 
gion while a coarser sampling may be used in the 
remaining part of the BZ (such non-uniform sam¬ 
pling was, however, not used in the present work). 
We stress that these facts apply only to isolated 2D 
semiconductors, which in practice means when a 
truncated Coulomb interaction is used. Only then 
ise(< 2 ) = \-{-0{q). If instead the full 1/r Coulomb 
interaction is used the calculations converge much 
faster but to a wrong value depending on the inter¬ 
layer distance.l^The dielectric function and Fig.[^ 


will be discussed in more depth in Sec. 3.6 



Figure 2: The static quasi-2D dielectric function 
of 2 H-M 0 S 2 along the F —)■ M direction. For com¬ 
parison the macroscopic dielectric function of bulk 
M 0 S 2 is also shown. The slope of the 2D dielectric 
function is indicated by a dashed line. 


Finally, we discuss the convergence of the 
GqWo energies with respect to the number of plane 
waves, Ng, used to represent the screened inter¬ 
action and self-energy. It has previously been 
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found that the GW corrections for bulk semicon¬ 
ductors and insulators follow a I/Nq scalingl^ 
which makes it possible to extrapolate the QP en¬ 
ergies to the infinite basis set limit. From our cal¬ 
culations with varying cut-off energy from 150eV 
and in some cases up to 500 eV we observe that: (i) 
The extrapolation procedure is essential and can 
correct QP energies obtained with 150eV cut-off 
by up to 0.5 eV. (ii) The slope of the extrapo¬ 
lation curve can be different for different states 
(bands and k-points), but generally shows a de¬ 
crease as function of Ng- (iii) The band gap tends 
to converge faster than the absolute band ener¬ 
gies. In Ref. |53]it was also shown that the lack 
of norm conservation of the PAW potentials can 
affect the convergence of the GW energies as Nq 
is increased. The effect is larger for more lo¬ 
calised states, particularly the states, where the 
violation of norm conservation can be significant. 
While it is possible to construct norm conserving 
PAW potentials, we have not pursued this in the 
present work. 


Hcut (eV) 



1/A/a 


Figure 3: GqWo quasi-particle energy of the va¬ 
lence band at the F point of monolayer 2 H-M 0 S 2 
as a function of 1 /Nq, where Nq is the number of 
plane waves. The different lines correspond to dif¬ 
ferent k-point samplings (/4 x x 1). The dashed 
lines shows the extrapolation to infinite plane wave 
cut-off. 

Performing the extrapolation to infinite cut-off 
for 30 X 30 k-points is computationally demand¬ 
ing. Fortunately, we have found that the extrapola¬ 
tion is rather insensitive to the k-point mesh. This 
is shown in Figure for the case of 2 H-M 0 S 2 : 


Changing the k-point mesh simply shifts the en¬ 
tire Ao-curves but do not affect their form. To ob¬ 
tain results converged with respect to both k-points 
and plane waves we have therefore performed the 
AG-extrapolation for a coarse k-point sampling of 
12 X 12 and corrected the band energies by the 
difference between a 12 x 12 and 30 x 30 cal¬ 
culation at 150eV plane wave cut-off. In doing 
this we have to interpolate band energies from the 
coarse to the fine k-point grid. The GW correc¬ 
tion to an energy, £k,„ on the fine grid is obtained 
as a weighted average of the corrections obtained 
at the nearest points of the coarse k-point grid, 
with weights determined by the overlap of 
the LDA state densities, 

where pkn(r) = |v/k«(r)p. Since the GqWo shift 
depends crucially on the shape of the wave func¬ 
tion, this approach is essential, in particular when 
interpolating the GqWq corrections close to band 
crossings. 


3 Results 

In this section we present the main results of our 
electronic structure calculations. To limit the pre¬ 
sentation, we have chosen to focus on the trends 
in electronic structure observed across the inves¬ 
tigated materials rather than giving in-depth anal¬ 
ysis of particular materials. However, as all the 
data is available the database it is straightforward 
for the interested reader to obtain the entire set of 
computed data. 

3.1 Stability 

The heat of formation of the 216 monolayer TMDs 
and TMOs in the relaxed structure was calculated 
from 


AA(MX2) = Ftot(MX2) - Fref (M) - 2Fref (X), 

( 1 ) 

where £'tot(MX 2 ) is the PBE total energy of the 
monolayer while Fref (M) and Aj-gf (X) are the refer¬ 
ence energies of the metal and chalcogen/oxygen 
(X), respectively. For the latter we use the fitted 
elemental phase reference energies (FERE).I^ 
Figure]^ shows the calculated heat of formations 
for all 216 compounds. Most of the materials have 
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negative heat of formation; in faet by requiring a 
heat of formation below 0.1 eV (to allow for uneer- 
tainties in the ealeulation methods) we obtain 171 
stable eompounds. The heat of formation of the 
stable semieonduetors together with relaxed lattiee 
eonstants (a), distanee between outermost ehaleo- 
gen/oxygen atoms (h) (see Figure [^, and the final 
magnetie moments are given in Table In gen¬ 
eral, the oxides have the highest stability followed 
by the sulphides, selenides, and tellurides in that 
order. Furthermore, the stability deereases as the 
metal ion goes through the transition metal series. 
For eomparison with previous studies we note that 
the 52 monolayer MX 2 eompounds found to be 
stable based on the LDA ealeulations of Ref. |36] 
form a subset of the stable materials identified in 
the present work. 


2H-MX2 
Sc 
Ti 
V 
Cr 
Mn 
Fe 
Co 
Ni 
Cu 
Ge 
Zr 
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Ru 
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Pb 

O S Se Te 
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Figure 4: Caleulated heat of formation for all 
monolayers in the 2H and IT phases. In general, 
the oxides have the highest stability followed by 
the sulphides, selenides, and tellurides in that or¬ 
der. Furthermore, the stability deereases as the 
metal ion goes through the transition metal series. 


While the heat of formation is a natural de- 
seriptor for whether the material will be possible 
to synthesize we stress that meehanieal instabil¬ 
ities or eompeting phases of lower energy have 
not been taken into aeeount. While it is possible 
to aeeount for both effeets, e.g. by earrying out 
moleeular dynamies simulations^^ and ineluding 
a larger pool of referenee systems!^ we have not 
pursued this further. Lebeque et al^ searehed the 
Inorganie Crystal Strueture Database (ICSD) for 
known layered bulk materials. They identified 46 
TMDs, but did not speeify whether the bulk ma¬ 
terials was known in the IT or 2H phase. Assum¬ 
ing the IT and 2H phases to be equally stable (it 
follows from Figure]^ that this is a reasonable as¬ 
sumption) 76 out of the 171 materials with nega¬ 
tive heat of formation, are already known as lay¬ 
ered bulk materials. These materials are marked 
by an asterix in Table The faet that Lebeque 
et al. did not identify the layered bulk form of 
any of the TMOs investigated here indieates that 
these struetures eould be very ehallenging to syn¬ 
thesise; presumable due to the existenee of non- 
layered bulk phases of higher stability. On the 
other hand, the monolayer TMOs might be meta¬ 
stable or eould be stabilised by interaetion with a 
substrate. Eneapsulation of the monolayers, as re- 
eently demonstrated for M 0 S 2 in hexagonal boron- 
nitride,!^ eould be a way to prevent the material 
from reaeting with other ehemieal speeies. 


3.2 Band gaps 


For all the stable and non-magnetie semieondue¬ 
tors we have performed GqWo ealeulations fol¬ 
lowing the proeedure deseribed in See. |2.2[ The 


GoWo eorreetions have been evaluated for the 10 
bands elosest to the Fermi energy and spin-orbit 
eoupling has been ineluded non-pertubatively. 
The spin-orbit splittings of the valenee/eonduetion 
bands for materials where these are non-vanishing 
are reported in Table Unless otherwise stated, 
all the results presented in this work inelude spin- 
orbit interaetions. As an example we show the 
GqWo and LDA band struetures of 2H-WSe2 in 
Figure For this partieular material we find a di- 
reet GqWo gap of 2.08 eV and a 0.45 eV splitting 
of the valenee band at the X-point. 

We refrain from providing a detailed eompari- 
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2 H-WSe2 



Figure 5: Band structure of 2H-WSe2 using LDA 
(black) and GqWo (blue) and LDA projected den¬ 
sity of states. Note the spin-orbit coupling gives 
rise to a splitting of the bands at various regions of 
the Brillouin zone. The blue line connecting the 
GqWo points is obtained from a cubic spline inter¬ 
polation. 


son with previous literature values for the QP band 
gaps. For M 0 S 2 such a comparison was made in 
Ref. |43l However, the fact that most previous GW 
calculations for 2D TMDs have used rather coarse 
k-point grids and have not employed a truncated 
Coulomb interaction (both of which have a sig¬ 
nificant effect on the calculated gap®l), the im¬ 
portance of spin-orbit interaction which is not al¬ 
ways included, as well as the sensitivity of the 
gap to the in-plane lattice constant (see Sec. |3.4[ ), 
makes a general comparison difficult. We believe 
that the QP band structures of the present work 
are the most carefully converged GqWo calcula¬ 
tions reported for 2D TMOs and TMDs to date. 
There are only few experimental reports on QP 
band gaps in freestanding TMD monolayers. For 
2 H-M 0 S 2 our GoWo gap of 2.48 eV agrees well 
with the 2.5 eV reported in Ref. Obased on pho¬ 
tocurrent spectroscopy on suspended M 0 S 2 . We 
note that the slightly smaller band gap obtained 
here for M 0 S 2 (2.48 eV) compared to our previous 
workl^(2.65eV) is mainly due to the inclusion of 
spin-orbit interaction in the present work. Alterna¬ 
tively, the QP gap can be inferred from optical ab¬ 
sorption or photo luminescence spectra which are 
experimentally simpler to obtain. However, this 
requires knowledge of the exciton binding energy 


which in turn depends on the screening from the 
substrate.^ 

In the last two columns of Table |3] we show 
the calculated direct and indirect band gaps of the 
51 stable 2D semiconductors. It is well known 
from photoluminescence spectroscopy that the 2H 
phase of monolayer M 0 S 2 , MoSe 2 , MoTe 2 , WS 2 , 
and WSe 2 have direct band gaps. This is repro¬ 
duced by our GqWo calculations. However we find 
that the indirect gap in MoTe 2 is about the same 
size as the direct gap. The only other materials we 
find to have a direct band gap are CrS 2 , CrSe 2 , and 
CrTe 2 with gaps of 1.54eV, 1.21 eV and O.VVeV, 
respectively. All other compounds have indirect 
gaps in the range 0.5-7.0 eV. 

In addition to LDA and GqWo we have cal¬ 
culated the band gaps using the GLLB-SC func¬ 
tional of Kuisma et al^ The GLLB-SC is an or¬ 
bital dependent exact exchange-based functional 
which, in addition to the Kohn-Sham band gap, 
provides an estimate of the derivative discontinu¬ 
ity. The GLLB-SC has previously been shown to 
yield good results for the band gap of bulk semi- 
conductorsbut to our knowledge it has not 
been previously applied to 2D materials. In Fig¬ 
ure we compare the GqWo band gaps with the 
PBE and GLLB-SC gaps. We first note that the 
GqWo band gaps range from 0.5 eV to almost 
8eV with the majority of the materials lying in 
the range leV to 3eV. We note that the size 
of the band gaps are directly correlated with the 
heat of formation of the materials with the ox¬ 
ides having the largest band gaps followed by the 
sulphides, selenides, and tellurides in that order. 
As expected, the LDA gaps are significantly lower 
than those obtained from GqWo which is consis¬ 
tent with the situation known from bulk materials 
and molecules. In contrast, except for a few out¬ 
liers, the band gaps obtained with the GLLB-SC 
functional lie very close to the GqWo values with 
a mean absolute error of 0.4 eV. This is consis¬ 
tent with the results obtained for both bulk and 
molecular systems^ and supports the use of the 
GLLB-SC functional as viable alternative to GW 
in large-scale studies where one would benefit of 
its low computational requirements that are similar 
to LDA. 
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Go Wo band gap (eV) 

Figure 6: Computed band gaps of 51 monolayer 
TMDs and TMOs. The GqWo band gaps are eom- 
pared to the band gaps obtained from DFT with 
the PBE and GLLB-SC xe-funetionals. The lat¬ 
ter ineludes the derivative diseontinuity of the xc- 
potential. Note the logarithmic scale. 

3.3 Absolute band positions 

For many applications not only the distance be¬ 
tween the occupied and unoccupied bands, i.e. the 
band gap, is of interest, but also the absolute posi¬ 
tion of the band edges relative to vacuum. We have 
calculated these by referring the band energies to 
the asymptotic value of the Hartree potential in the 
vacuum region between the layers. For bulk mate¬ 
rials this is a difficult task as it requires the use 
of thick slabs to represent both the bulk interior 
and its surface. Moreover the Hartree potential 
depends on the surface dipoles (on both sides of 
the slab) which makes the problem highly surface 
dependent and complicates the comparison with 
experiments. These problems are obviously not 
present for the the monolayers studied here mak¬ 
ing them ideal as benchmarking systems for the 
band alignment problem. 

In Figure |7] the positions of the valence band 
maximum (VBM) and conduction band minimum 
(CBM) relative to the vacuum level are shown for 
the different oxides and chalcogenides at both the 
LDA and GqWo level. As a significant part of 
the GLLB-SC band gap comes from the deriva¬ 
tive discontinuity which applies to the fundamen¬ 
tal gap rather than the individual band energies, the 
GLLB-SC cannot be used to obtain the absolute 
band edge positions. For all materials, the effect 


of the Go Wo correction is to shift the conduction 
band up and the valence band down with respect 
to the LDA values. In fact, the corrections of the 
VBM and CBM are rather symmetric meaning that 
the band gap centre is largely unaffected by the 
GoWo correction (see below). 

It has been suggested that 2D semiconductors 
could be used for photo-catalytic water splitting. 
This is mainly motivated by their excellent light 
absorption, large specific surface area, and read¬ 
ily tuneable electronic properties The equi¬ 
librium potentials for the hydrogen- and oxygen 
evolution reactions at pH 7 are indicated by dashed 
green lines in Figure Materials with CBM 
above the standard hydrogen electrode (SHE) at 
—4.03eV relative to vacuum (at pH 7) , could in 
principle be used to evolve hydrogen at the cath¬ 
ode of a photo-catalytic water splitting device.l^ 
Likewise materials with VBM below the oxygen 
evolution potential (1.23 eV below the SHE) could 
in principle be used a photo anode in the wa¬ 
ter splitting reaction. In practice, the CBM/VBM 
should lie a few tenths of an eV above/below the 
redox potentials to account for the intrinsic en¬ 
ergy barriers of the water-splitting reactions re¬ 
actions.® As can be seen a number of the TMD 
monolayers qualify as potential water splitting 
photo-electrodes based on their energy level posi¬ 
tions. A very citical issue, however, is the stability 
of the materials under the highly oxidizing reac¬ 
tion conditions. A possible solution to this prob¬ 
lem could be to protect the photo-electrode from 
direct contact with the water by a transparent and 
highly stable thin film, which in practice means an 
oxide material. 

A simple empirical relation between the band 
gap center of a semiconductor and the electroneg¬ 
ativities of the constituent atoms has been sug¬ 
gested:® 

£ce„ter=-[z(M)z(X)2]l/3, (2) 

where 2 ^(M) and xO^) the electronegativity 
of the metal and oxygen/chalcogen on the Mul- 
liken scale, respectively. In Figure we compare 
the band gap centers obtained from GqWo with 
those obtained from LDA and calculated with 
Eq. ??, where experimentally obtained values of 
the electronegativities® have been used. The 
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Figure 7: Position of the valence band maximum and conduction band minimum relative to the vacuum 
level (set to zero) for LDA and GqWo . In both cases spin-orbit splitting of the bands has been taken into 
account. The hydrogen and oxygen evolution potentials at pH 7 are shown by green dashed lines. 


band gap centers from LDA and GqWo agree quite 
well showing a mean absolute deviation from the 
GqWo values of only 0.2eV. While it is known 
that the Kohn-Sham band gap center is formally 
exact within DFT,l^it is somewhat surprising that 
the LDA performs that well. While the empirical 
formula is able to describe the qualitative trends 
of the gap centres the quantitative values devi¬ 
ate significantly from the ab-initio results, with a 
mean absolute difference from the GqWo result 
of 0.9eV and a mean relative deviation of 14%. 
We ascribe a large part of this deviation to origi¬ 
nate from dipole fields formed due to the positively 
charged metal ions and negatively charge chalco- 
gens/oxygens which will increase potential outside 
the monolayer and thereby down shift the bands - 
an effect not accounted for by the empirical for¬ 
mula. Since the size of the dipoles is determined 
by the amount of charge transfer, the deviation be¬ 
tween Eq. (??) and the ab-initio results is expected 
to correlate with the difference in electronegativ¬ 
ity between the metal atom and chalcogen/oxygen 
atoms. From the inset of Figure we see that 
this indeed is the case: For materials with larger 
difference in electronegativity between the atomic 
species (A;^) the band gap center given by Eq. (??) 
generally deviates more from the GqWo results. 

While it is important to establish the intrinsic 
properties of the 2D materials in their isolated 
form, practical applications as well as most ex¬ 
perimental setups, involve heterostructures where 



Figure 8: Comparison of the absolute band gap 
centers (relative to vacuum) obtained from GqWo , 
LDA, and the empirical formula Eq. (??). In¬ 
set: The difference in the band gap centers from 
GoWo and Eq. (??), AEeen = com¬ 

pared to the difference in the electronegativities 
of the metal and oxygen/chalcogen atom. Ax = 
X{M)-x{X). 

the 2D materials are stacked into van der Waals 
heterostructure or simply lie on a substrate. In 
such systems the alignment of the bands at the het¬ 
erostructure interfaces become crucial. Due to the 
weak interaction between 2D semiconductors it is 
reasonable to expect that the band alignment at the 
interface between two different 2D can be obtained 
by aligning the band edges of the isolated systems 
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relative to a eommon vaeuum level. This is equiv¬ 
alent to disregarding effeets of band hybridiza¬ 
tion and the formation of interfaee dipoles due to 
eharge redistribution. Verifying this assumption 
from first-prineiples ealeulations is, however, dif- 
fieult due to the lattiee mismateh between different 
2D materials. 

To provide an overview of the band edge posi¬ 


tions of the 51 monolayers, we show in Fig. 9a 


the CBM plotted against the VBM obtained from 
GoWo . To illustrate the use of sueh a diagram we 
have highlighted 2 H-M 0 S 2 and indieated regions 
eorresponding to different band alignments with 
M 0 S 2 . The possible band alignments are: Strad¬ 
dling gap (type I), staggered gap (type II), and bro¬ 
ken gap (type III). For many applieations, e.g. tan¬ 
dem photovoltaie deviees or ereation of long lived 
indireet exeitons, a type II band alignment is pre¬ 
ferred. We have highlighted a few materials that 
are expeeted to form type II band alignment with 
M 0 S 2 . The detailed band alignments for these ma¬ 


terials are shown in Fig. 9b 


3.4 Strain effects on the band struc¬ 
ture 

In the present work we have eonsidered monolay¬ 
ers in their PBE relaxed geometry. Sinee PBE er¬ 
rors on lattiee eonstants typieally are around 1 to 
2% we have investigated how this would affeet the 
LDA band strueture. From Figure [T^ we see that a 
ehange of the lattiee parameter within the eonsid¬ 
ered range ean produee quite drastie ehanges in the 
band gap. For example, in the ease of 2 H-M 0 S 2 
a ehange in the lattiee eonstant from the PBE 
value (3.18 A) to the experimental value (3.16A) 
ehanges the band gap by around 0.1 eV. From Ta¬ 
ble]^ where the band gaps are given as funetion of 
lattiee eonstant, we furthermore see that the LDA 
gap ehanges from indireet to direet under 1 % eom- 
pressive strain. A few other direet gap materials 
are seen to develop an indireet gap when strained. 
Thus we eonelude that both the size and nature of 
the band gap of the monolayers ean depend deli- 
eately on the lattiee eonstant. 

To understand the different behaviour of the 
band gap upon strain, we have analysed the pro- 
jeeted density of states (see supplementary mate¬ 
rial). We find that the materials ean be roughly 


divided into two elasses aeeording to the nature 
of the wave funetions around the band gap. For 
the materials with group 6 metals (Cr, Mo and 
W), the valenee and eonduetion band states are 
bonding/anti-bonding eombinations of the metal 
J-states and oxygen/ehaleogen p-states and in their 
equilibrium lattiee eonstant they have direet band 
gaps. For these materials we find that inereas- 
ing the lattiee eonstant, inereases the M-X bind¬ 
ing distanee whieh weakens the hybridization and 
reduees the bonding/anti-bonding gap. The other 
elass is TMDs with metals from group 4, 10 and 
14 (Ti, Zr, Hf, Ni, Pd, Pt, Ge, Sn, Pb). For these 
materials, the valenee band states have primarily 
ehaleogen /)-eharaeter while the eonduetion band 
is either metal-J (group 4), ehaleogen-/? (group 
10) or metal-5 and ehaleogen-/? (group 14). In 
these eases the gap size is eontrolled by the width 
of the eonduetion band and the ehaleogen valenee 
band; applieation of a tensile strain will eause the 
states to beeome more loealized narrowing eaeh 
of the bands and thereby opening the gap. As a 
eonsequenee of the deeoupled bands, these also all 
have indireet band gaps. 

3.5 Effective masses 

From the GqWo band struetures we have extraeted 
the effeetive eleetron and hole masses by fitting 
a paraboloid to the energies of the 19 nearest 
k-points around the eonduetion band minimum 
(CBM) and valenee band maximum (VBM) ae¬ 
eording to 

E = — k^Ak, (3) 

ZIHq 

where k = (kx, ky) is the in plane k-point measured 
from the band extremum. The eigenvalues of the 
matrix A yield the inverse effeetive masses in the 
direetion of smallest and largest eurvature. If the 
CBM or VBM is loeated at one of the high symme¬ 
try points of the BZ (the F or X point) the effeetive 
masses will naturally be isotropie. However, for 
band extrema loeated at other points this is gener¬ 
ally not the ease. 

In Figure [TTa] and |b] we show the effeetive elee¬ 
tron and hole masses along the two natural di- 
reetions. Points falling on the diagonal line eor- 
respond to isotropie band masses. The effeetive 
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Figure 9: Band alignment diagram, (a) Conduction band minimum ^cbm plotted against the valence band 
maximum F'vbm for the 51 monolayers. The band edges relative to vacuum are obtained from GqWo . 
As an example we have highlighted 2 H-M 0 S 2 (orange dot) and indicated the regions corresponding to the 
different types of band alignment: Straddling gap (type I), staggered gap (type II), and broken gap (type 
III). A few selected materials that will form type-II heterostructures with M 0 S 2 have been highlighted in 
green, (b) Absolute band edge positions and band gaps of 2 H-M 0 S 2 and the selected materials highlighted 
in (a). 


masses are also listed in Table [51 We note that 
the effective electron masses lie in the range 0.1 
to lOme with roughly an equal number being light 
(m* < me) and heavy (m* > m^). The same approx¬ 
imately applies to the hole masses, although they 
seem to be generally heavier than the electrons. 
In accordance with the discussion in the previous 
section, we see that only the materials with direct 
gaps (group 6 metals) have both isotropic electron 
and hole masses. For other materials the masses 
can be quite anisotropic and we would also expect 
the masses to depend sensitively on the lattice con¬ 
stant. 

In order to estimate exciton binding energies 


(see Sec. |3.7| ) we also evaluate the effective ex¬ 
citon masses defined as 


1 . — 1 * — 1 I * — 1 

Mex ="*e +'«h ■ 


(4) 


We distinguish between two kinds of excitons: di¬ 
rect excitons that possess zero momentum and in¬ 
direct gap excitons that have a finite momentum 
corresponding to the distance in k-space between 


the VBM and CBM. In Fig. Figure jllcj we plot 
the effective exciton masses along the two natural 
directions (we show both the direct and indirect 


exciton mass whether the material has direct or in¬ 
direct gap). 


3.6 Dielectric function 


The dielectric function is one of the most impor¬ 
tant material response functions. It relates the 
strength of an externally applied field to the to¬ 
tal (screened) field in the material. In particular, 
it determines the strength of the electron-electron 
interaction and is a key ingredient in calculations 
of electronic states such as QP band structures and 
excitons. 

For many purposes it is not necessary to know 
the precise spatial variation of the induced poten¬ 
tials but only its average value over a unit cell. The 
relation between the external potential and the av¬ 
eraged total potential is described by the macro¬ 
scopic dielectric function which can be obtained 
from the microscopic dielectric function according 


to. 


1 

£M(q, co) 




(5) 


Here egg,(q, m) is the plane wave representa¬ 
tion of the inverse microscopic dielectric func- 
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Figure 10: Change in the LDA band gap when the in-plane lattice constant is varied between —2% and 
2 %. 
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Figure 11: Effective electron (a), hole (b) and exciton (c) masses (in units of me) along the two princi¬ 
pal directions obtained. The masses are calculated from the GqWo band structures including spin-orbit 
interaction. Points on the dashed lines correspond to isotropic masses. 


tion which is a standard output of many electronic 
structure codes. For bulk semiconductors one usu¬ 
ally refers to the ^ = 0 and to = 0 limit of Cm as 
the dielectric constant. 

In the case of a 2D material Eq. (??) must be 
generalized as there is no natural unit cell over 
which to perform the average of the total field. If 
one restricts the averaging region to a slab of width 
d containing the 2D material one arrives at the fol¬ 
lowing expression:®! 





(qii)- (6) 


We note that due to the averaging procedure, 
takes the finite thickness of the material into ac¬ 
count. We therefore refer to it as a quasi-2D di¬ 


electric function to distinguish it from a mathe¬ 
matically strict 2D quantity where the third dimen¬ 
sion has been integrated out. As discussed in Ref. 
I43l e^(q[| = 0) = 1 which implies that long wave 
length perturbations are not screened by the 2D 
material at all. In particular, there is no direct ana¬ 
logue of the dielectric constant in 2D; any realis¬ 
tic model for screening in 2D materials must be 
( 2 '-dependent. 

We have calculated e^((?||) along the F —M 
and r .ST directions for the 51 stable 2D semi¬ 
conductors. We have found that this quantity is al¬ 
most isotropic within the plane of the monolayer. 
The thickness of the averaging region has been set 
to d = 2h, where h is the thickness of the layer, but 
as show in Ref. @3| the dielectric function is not 
very sensitive to this value - in particular for the 
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most important regime of < 2 ' < \/d. 

As an example Figure]^ shows the statie maero- 
seopie dieleetrie funetion of 2 H-M 0 S 2 . The linear 
inerease for small q followed by a maximum and 
then a monotonie deerease towards 1 in the large 
q limit is eharaeteristie for all 2D semieonduetors. 
For eomparison we also show the dieleetrie fune¬ 
tion of bulk M 0 S 2 for the same in-plane q veetors. 
To tie up with the diseussion in See. |2.2| we note 
that it is the strong ^-dependenee of for small 
q that is responsible for the very slow k-point eon- 
vergenee of the GW ealeulations. 

To illustrate the variation in the dieleetrie proper¬ 
ties of the monolayers we show the slope of {q) 
at (? = 0 in Figure Not surprisingly, the varia¬ 
tion eorrelates well with the size of the eleetronie 
band gaps, also shown in the figure: Large band 
gap materials have smaller dieleetrie funetion and 
viee versa. The slopes of the maeroseopie dielee¬ 
trie funetion are listed in Tabled 


3.7 Excitons 


One of the most eharaeteristie features of atomi- 
eally thin 2D semi eondueto rs is the large binding 
energy of exeitons.^^^^^^^ The reason for this is 
the redueed sereening due to the lower dimension 
whieh yields a stronger attraetion between elee- 
trons and holes (see diseussion in previous see- 
tion). The eonventional method for ealeulating ex- 
eiton binding energies from first prineiples is the 
Bethe-Salpeter equation (BSE). The BSE is eom- 
putationally highly demanding and not suited for 
large-seale studies like the present. Instead we 
use a reeently developed 2D Mott-Wannier model 
for exeitons that only needs the exeiton effeetive 
mass and the quasi-2D dieleetrie funetion as in¬ 
put. In real spaee the model takes the form of a 2D 
Sehrodinger equation. 




V/(r) =EbV^(r), 


(7) 


where /iex is the effeetive exeiton mass and W{r) 
is the l/r Coulomb interaetion between the elee- 
tron and the hole sereened by the non-loeal e^. 
The model has been benehmarked against full BSE 
ealeulations for 2 H-M 0 S 2 and 2 H-WS 2 and the re¬ 
sults were found to deviate by less than 0.1 eV. 


The four basie assumptions behind the Mott- 
Wannier exeiton model are: (i) Isotropie exeiton 
masses, (ii) parabolie band struetures elose to the 
fundamental gap, (iii) the exeiton is well deseribed 
by transitions between the valenee and eondue- 
tion band only, and (iv) the valenee and eondue- 
tion band wave funetions are uniformly distributed 
over the layer, i.e. their profile along z ean be ap¬ 
proximated by a step funetion. While the dieleetrie 
funetions were found to be very nearly isotropie 
for all materials, this is not the ease for the exei¬ 
ton masses, see Eigure [TT] (e). While it is possi¬ 
ble to modify the model to allow for anisotropie 
masses we here limit ourselves to the materials 
with isotropie exeiton masses. The exeiton bind¬ 
ing energies obtained from the model are shown 
as the dark region on the top of the bars in Eig¬ 
ure 13 see also Table The total height of the 
bar represents the GoWq ealeulated QP gap. Eor 
direet (indireet) band gap materials we have used 
the direet (indireet) exeiton mass in the model. 

In aeeordanee with earlier experimental and the- 
oretieal studies we find strong exeiton binding en¬ 
ergies on the order of 20-30% of the band gap. In 
general materials with larger QP band gaps have 
more strongly bound exeitons. This follows from 
the eorrelation between the size of the band gap 
and the dieleetrie funetion in Eigure Earger 
band gap implies a smaller dieleetrie funetion and 
thus a stronger eleetron-hole interaetion. In Ta¬ 
ble we eompare our ealeulated exeiton binding 
energies with optieal data from experiments. We 
find good agreement for MoSj, MoSe 2 and WSe 2 
while the agreement is less satisfaetory for MoTe 2 
and WS 2 . It should be noted, however, that the ex¬ 
perimental exeiton binding energy for MoTe 2 was 
obtained as the differenee between the ealeulated 
GqWo band gap and the position of the optieal 
photoluminiseenee peak. Thus inaeeuraeies in the 
GqWo band gap as well as substrate effeets on the 
measured photolumineseenee peak eould explain 
the disagreement. 


4 Conclusions 

We have presented a detailed eleetronie strueture 
study of 51 monolayer transition metal diehaleo- 
genides and -oxides. The 51 monolayers were 
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Figure 12: The slope of the static quasi-2D dielectric function, e^(^), evaluated dXq = 0. The materials 
are ordered according to their LDA direct band gap. 


Table 6: Exciton binding energies in eV calcu¬ 
lated from the Mott-Wannier model compared 
to experimental values. 


name Eh (model) Eh (exp.) 


2 H-M 0 S 2 0.47 

2H-MoSe2 0.42 

2H-MoTe2 0.36 

2 H-WS 2 0.48 

2H-WSe2 0.43 


O. 55 E] 

0.5^21 

ObISSb 

0.66,1^ 0.7 1 I 22 ! 
0.38, El 0 . 37 !^ 


“The exciton binding energy is obtained from subtracting 
the energy of the measured exciton photoluminiscence peak 
from our calculated Go Wo band gap. 


chosen out of an initial set of 216 compounds as 
those having a finite band gap and a negative heat 
of formation. The calculated properties include the 
LDA band structure for in-plane lattice constants 
in a range around the equilibrium structure, the 
quasiparticle band structure at the equilibrium lat¬ 
tice constant evaluated in the GqWo approximation 
and including spin-orbit coupling, the absolute po¬ 
sitions of the conduction and valence band edges 
relative to vacuum, the effective electron and hole 
masses, and the static ^-dependent dielectric func¬ 
tions. As an example we showed how the com¬ 
puted data, in this case the effective masses and di¬ 
electric functions, can be used to obtain the lowest 
exciton binding energies from a 2D Mott-Wannier 
model. 

Rather than providing a detailed account of the 
electronic structure of specific materials, we have 



Figure 13: GqWo band gaps (total bar height) and 
exciton binding energies (darker topmost part of 
the bar). The exciton binding energy was obtained 
from a quasi-2D Mott-Wannier model. Only ma¬ 
terials with isotropic exciton masses are shown. 
The green and blue bars refer to indirect and di¬ 
rect band gaps and excitons, respectively. 


chosen to focus on general trends and correlations 
in the electronic structure of the materials. How¬ 
ever, as all the computed data is available in an 
open database it is straightforward to retrieve and 
analyse specific materials data in greater detail. 
We are presently working to expand the database 
to include other 2D materials and properties. We 
strongly believe that such a database will be use¬ 
ful both for guiding experimental efforts in the the 
search for new 2D materials and as a platform for 
predicting properties of more complex materials 
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such as van der Waals heterostructures. 

5 Supporting information 

Figures showing the LDA and GqWo band strue- 
tures of all 51 stable non-magnetie semieondue- 
tors and projeeted density of states. This mate¬ 
rial is available free of eharge via the Internet at 

http://pubs.acs.org, 
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Table 1: Relaxed in-plane lattice constant, a, distance between chalcogen/oxygen atoms, h, forma¬ 
tion energies from PBE, and using tbe fitted elemental phase reference energies, and 

total magnetic moment, /i. A * denotes whether the material is found in bulk form according to Ref. 

m 


name 

a (A) 

/1(A) 

(eV) 

£:FERE (gY) 

H (Mb) 

2H-Cr02 

2.63 

2.34 

-1.64 

-1.99 

0.0 

2H-CrS2* 

3.05 

2.95 

-0.662 

-0.892 

0.0 

2H-CrSe2* 

3.21 

3.15 

-0.474 

-0.65 

0.0 

2H-CrTe2* 

3.47 

3.41 

-0.051 

-0.104 

0.0 

2H-Ge02 

2.81 

2.32 

-0.969 

-1.28 

0.0 

1T-Ge02 

2.9 

1.96 

-1.53 

-1.84 

0.0 

lT-GeS 2 

3.44 

2.8 

-0.222 

-0.416 

0.0 

2H-Hf02 

3.12 

2.34 

-2.71 

-3.08 

0.0 

1T-Hf02 

3.25 

1.95 

-3.27 

-3.64 

0.0 

2H-HfS2* 

3.54 

3.14 

-1.37 

-1.62 

0.0 

lT-HfS2* 

3.65 

2.9 

-1.59 

-1.83 

0.0 

2H-HfSe2* 

3.68 

3.36 

-1.17 

-1.36 

0.0 

IT-HfSez* 

3.77 

3.16 

-1.34 

-1.53 

0.0 

2H-HfTeo* 

3.91 

3.7 

-0.656 

-0.723 

0.0 

1T-Mn02 

2.89 

1.93 

-1.58 

-2 

3.0 

2 H-M 0 O 2 

2.82 

2.45 

-1.73 

-1.94 

0.0 

2 H-M 0 S 2 * 

3.18 

3.13 

-0.842 

-0.93 

0.0 

2H-MoSe2* 

3.32 

3.34 

-0.663 

-0.698 

0.0 

2H-MoTe2* 

3.55 

3.61 

-0.237 

-0.149 

0.0 

1T-Ni02 

2.84 

1.91 

-0.716 

-1.01 

0.0 

lT-NiS2 

3.35 

2.35 

-0.248 

-0.424 

0.0 

lT-NiSe 2 

3.54 

2.49 

-0.251 

-0.374 

0.0 

1T-Pb02 

3.39 

2.14 

-0.641 

-0.8 

0.0 

lT-PbS2 

3.85 

3.09 

0.069 

0.031 

0.0 

1T-Pd02 

3.09 

1.96 

-0.272 

-0.482 

0.0 

lT-PdS2* 

3.55 

2.49 

-0.125 

-0.214 

0.0 

IT-PdSez* 

3.73 

2.63 

-0.206 

-0.242 

0.0 

lT-PdTe 2 ‘ 

4.02 

2.76 

-0.177 

-0.09 

0.0 

1T-Pt02 

3.14 

1.9 

-0.405 

-0.612 

0.0 

lT-PtS2* 

3.57 

2.46 

-0.332 

-0.418 

0.0 

lT-PtSe 2 * 

3.75 

2.62 

-0.364 

-0.397 

0.0 

lT-PtTe 2 * 

4.02 

2.77 

-0.321 

-0.23 

0.0 

2H-Sc02 

3.22 

2.07 

-2.37 

-2.74 

1.0 

2H-ScS2 

3.79 

2.72 

-1.21 

-1.46 

1.0 

2H-ScSe2 

3.95 

2.94 

-1.1 

-1.29 

1.0 

2H-Sn02 

3.09 

2.46 

-0.225 

-0.54 

0.0 

1T-Sn02 

3.22 

2 

-1.01 

-1.33 

0.0 

2H-SnS2* 

3.61 

3.23 

-0.048 

-0.241 

0.0 

lT-SnS 2 * 

3.7 

2.96 

-0.333 

-0.527 

0.0 

lT-SnSe 2 * 

3.86 

3.19 

-0.285 

-0.425 

0.0 

2H-Ti02 

2.88 

2.26 

-1.83 

-2.02 

0.0 

1T-Ti02 

2.99 

1.94 

-2.91 

-3.1 

0.0 

2H-TiS2* 

3.34 

3.02 

-1.16 

-1.23 

0.0 

2H-TiSe2* 

3.49 

3.24 

-1 

-1.02 

0.0 

2H-TiTe2‘ 

3.74 

3.58 

-0.544 

-0.441 

0.0 

2H-VSe2* 

3.34 

3.2 

-0.699 

-0.956 

1.0 

2H-VTe2* 

3.6 

3.5 

-0.263 

-0.397 

1.0 

2 H-WO 2 

2.83 

2.48 

-1.74 

-1.85 

0.0 

2 H-WS 2 * 

3.19 

3.15 

-0.783 

-0.776 

0.0 

2H-WSe2* 

3.32 

3.36 

-0.547 

-0.487 

0.0 

2H-Zr02 

3.14 

2.33 

-2.65 

-2.96 

0.0 

1T-Zr02 

3.26 

1.93 

-3.18 

-3.49 

0.0 

2H-ZrS2* 

3.57 

3.14 

-1.37 

-1.55 

0.0 

lT-ZrS 2 * 

3.68 

2.9 

-1.55 

-1.47 

0.0 

2H-ZrSe2* 

3.7 

3.37 

-1.2 

-1.33 

0.0 

lT-ZrSe 2 * 

3.79 

3.16 

-1.34 

-1.47 

0.0 

2H-ZrTe2’ 

3.92 

3.73 

-0.739 

-0.746 

0.0 
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Table 2: Spin orbit induced splittings at the valence band maximum and conduction band minimum 
as found in the LDA and Go Wo band structure respectively. Materials with negligible spin-orbit 
coupling are not shown. The location of the band extremum in the BZ is indicated in parenthesis®. 
Note that this can be different in LDA and GoWo . 


name 

(LDA) 

(LDA) 

(GoWo) 

A£c^b°m (GoWo) 

2H-CrS2 

0.07 (K) 

0(K) 

0.07 (K) 

0.01 (K) 

2H-CrSe2 

0.09 (K) 

0.02 (K) 

0.1 (K) 

0.02 (K) 

2H-CrTe2 

0.12 (K) 

0.02 (K) 

0.13 (K) 

0.03 (K) 

2H-Hf02 

0(T) 

0.17 (T) 

0(T) 

0.15 (T) 

2H-HfS2 

0.03 (T) 

0.07 (T) 

0.02 (T) 

0.09 (X) 

2H-HfSe2 

0.13 (T) 

0.1 (T) 

0.12 (T) 

0.11 (X) 

2 H-Hfre 2 

0(r) 

0.15 (T) 

0.48 (T) 

0.18 (X) 

2H-M0S2 

0.15 (K) 

0(K) 

0.15 (K) 

0(K) 

2H-MoSe2 

0.19 (K) 

0.02 (K) 

0.19 (K) 

0.02 (K) 

2H-MoTe2 

0.23 (K) 

0.04 (K) 

0.25 (K) 

0.05 (X) 

2H-Ti02 

0(T) 

0.02 (T) 

0(X) 

0.02 (T) 

2H-TiS2 

0.02 (T) 

0(T) 

0(r) 

0(E) 

2H-TiTe2 

0(r) 

0(r) 

0.32 (T) 

0(E) 

2H-WO2 

0(r) 

0.02 (T) 

0(r) 

0(K) 

2H-WS2 

0.45 (K) 

0.04 (K) 

0.45 (K) 

0.02 (K) 

2H-WSe2 

0.49 (K) 

0.04 (K) 

0.49 (K) 

0.03 (K) 

2H-Zr02 

0(T) 

0.05 (T) 

0(T) 

0.05 (T) 

2H-ZrS2 

0.02 (T) 

0.02 (T) 

0.02 (T) 

0(E) 

2H-ZrSe2 

0.1 (T) 

0.03 (T) 

0.1 (T) 

0(E) 

2H-ZrTe2 

0(r) 

0(r) 

0.28 (T) 

0(E) 


^See Fig. 
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Table 3: Absolute band edge positions with respect to vacuum, their location in the Brillouin zone 
(in parenthesis)" and corresponding band gaps as obtained by LDA and GqWo . 


LDA GoWo@LDA 


name 

•^vbm 

(eV) 

^cbm 

(eV) 

£gap 

(eV) 

pdirect 

-^gap 

(eV) 

^vbm 

(eV) 

•^cbm 

(eV) 

^gap 

(eV) 

rdirect 

-^gap 

(eV) 

2H-Cr02 

-7.64 (P) 

- 7.21 (K) 

0.43 

1.57 (K) 

- 7.37 (P) 

- 5.73 (K) 

1.64 

2.45 (K) 

2H-CrS2 

-6.08 (K) 

-5.18 (K) 

0.90 

0.90 (K) 

-5.85 (K) 

-4.31 (K) 

1.54 

1.54 (K) 

2H-CrSe2 

-5.50 (K) 

-4.80 (K) 

0.70 

0.70 (K) 

-5.22 (K) 

-4.02 (K) 

1.21 

1.21 (K) 

2H-CrTe2 

-5.05 (K) 

-4.59 (K) 

0.45 

0.45 (K) 

-4.75 (K) 

-3.98 (K) 

0.77 

0.77 (K) 

2H-Ge02 

-8.99 (K) 

-7.62 (P) 

1.37 

1.95 (P) 

-10.65 (X) 

-6.40 (P) 

4.24 

4.62 (P) 

1T-Ge02 

-9.20 (T) 

-5.65 (P) 

3.55 

4.04 (P) 

-11.07(T) 

-4.00 (P) 

7.07 

7.55 (P) 

lT-GeS 2 

-6.63 (E) 

-6.07 (M) 

0.57 

1.00 (M) 

-7.57 (E) 

-5.38 (M) 

2.19 

2.64 (M) 

2H-Hf02 

-8.51 (T) 

-6.70 (T) 

1.80 

1.91 (T) 

-9.91 (T) 

-5.41 (T) 

4.50 

4.60 (T) 

1T-Hf02 

-8.24 (T) 

-3.61 (M) 

4.63 

4.85 (E) 

-9.89 (X) 

-1.91 (M) 

7.98 

8.21 (E) 

2H-HfS2 

-7.05 (T) 

-6.12 (X) 

0.93 

1.20 (X) 

-8.14 (T) 

-5.52 (X) 

2.63 

2.93 (X) 

lT-HfS2 

-6.48 (P) 

-5.42 (M) 

1.06 

1.77 (P) 

-7.62 (P) 

-4.63 (M) 

2.98 

3.97 (P) 

2H-HfSe2 

-6.47 (T) 

-5.77 (X) 

0.70 

1.02 (X) 

-7.38 (T) 

-5.29 (X) 

2.09 

2.49 (X) 

lT-HfSe 2 

-5.57 (P) 

-5.26 (M) 

0.30 

1.08 (P) 

-6.53 (P) 

-4.58 (M) 

1.96 

2.95 (P) 

2H-HfTe2 

-5.46 (P) 

-5.39 (X) 

0.06 

0.52 (E) 

-6.06 (T) 

-5.12 (X) 

0.94 

1.62 (X) 

2 H-M 0 O 2 

-6.99 (P) 

-6.09 (K) 

0.91 

1.66 (P) 

-7.37 (P) 

-5.17 (K) 

2.20 

2.94 (P) 

2 H-M 0 S 2 

-6.13 (K) 

-4.55 (K) 

1.58 

1.58 (K) 

-6.32 (K) 

-3.84 (K) 

2.48 

2.48 (K) 

2H-MoSe2 

-5.50 (K) 

-4.18 (K) 

1.32 

1.32 (K) 

-5.63 (K) 

-3.46 (K) 

2.18 

2.18 (K) 

2H-MoTe2 

-5.04 (K) 

-4.11 (K) 

0.93 

0.93 (K) 

-5.11 (K) 

-3.40 (X) 

1.71 

1.72 (K) 

1T-Ni02 

-8.38 (X) 

-7.22 (E) 

1.17 

1.37 (E) 

-8.38 (T) 

-6.24 (E) 

2.15 

2.31 (E) 

lT-NiS2 

-5.97 (P) 

-5.46 (E) 

0.51 

0.89 (T’) 

-6.61 (E) 

-4.24 (E) 

2.38 

2.76 (X) 

IT-NiSez 

-5.10 (P) 

-5.10(E) 

0.00 

0.56 (T’) 

-5.70 (E) 

-3.91 (E) 

1.79 

2.24 (X) 

1T-Pb02 

-8.47 (X) 

-7.15 (P) 

1.32 

1.58 (P) 

-9.50 (X) 

-6.47 (P) 

3.03 

3.26 (P) 

lT-PbS2 

-6.93 (X) 

-6.29 (M) 

0.63 

0.81 (M) 

-7.67 (E) 

-5.95 (M) 

1.72 

1.91 (M) 

1T-Pd02 

-7.82 (X) 

-6.52 (E) 

1.30 

1.71 (E) 

-8.20 (E) 

-5.36 (E) 

2.84 

3.24 (E) 

lT-PdS2 

-6.48 (P) 

-5.37 (E) 

1.11 

1.30(E) 

-7.19 (T) 

-4.70 (E) 

2.48 

2.65 (X) 

IT-PdSej 

-5.56 (P) 

-5.08 (E) 

0.48 

0.87 (T’) 

-6.25 (E) 

-4.46 (E) 

1.79 

2.10(E) 

lT-PdTe 2 

-4.46 (P) 

-4.62 (E) 

0.00 

0.40 (T’) 

-5.20 (P) 

-4.18(E) 

1.02 

1.39 (T’) 

1T-Pt02 

-7.21 (X) 

-5.61 (E) 

1.60 

2.00 (E) 

-7.99 (E) 

-4.41 (E) 

3.59 

4.00 (E) 

lT-PtS2 

-6.44 (T) 

-4.84 (E) 

1.61 

1.69(E) 

-7.16(E) 

-4.21 (E) 

2.95 

3.14 (T’) 

lT-PtSe 2 

-5.69 (P) 

-4.62 (E) 

1.07 

1.29 (E) 

-6.52 (T) 

-4.04 (E) 

2.48 

2.67 (E) 

lT-PtTe 2 

-4.52 (P) 

-4.29 (E) 

0.23 

0.75 (T’) 

-5.44 (P) 

-3.74 (E) 

1.69 

2.03 (T’) 

2H-Sn02 

-8.78 (K) 

-8.21 (P) 

0.56 

1.26 (P) 

-10.16 (X) 

-7.50 (P) 

2.66 

3.31 (P) 

1T-Sn02 

-8.64 (X) 

-6.10 (P) 

2.54 

3.13 (E) 

-10.27 (X) 

-4.89 (P) 

5.38 

5.93 (E) 

2H-SnS2 

-6.54 (P) 

-5.95 (M) 

0.59 

0.91 (P) 

-7.54 (P) 

-5.61 (M) 

1.93 

2.14 (P) 

lT-SnS 2 

-6.98 (X) 

-5.58 (M) 

1.40 

1.65 (M) 

-7.98 (E) 

-4.91 (M) 

3.07 

3.33 (M) 

lT-SnSe 2 

-6.19 (P) 

-5.58 (M) 

0.62 

0.96 (M) 

-6.96 (E) 

-5.05 (M) 

1.91 

2.25 (M) 

2H-Ti02 

-8.88 (T) 

-7.78 (T) 

l.IO 

1.25 (X) 

-9.97 (X) 

-6.25 (T) 

3.72 

3.83 (X) 

1T-Ti02 

-8.67 (X) 

-6.02 (P) 

2.65 

2.80 (P) 

-9.80 (X) 

-4.07 (E) 

5.74 

5.97 (E) 

2H-TiS2 

-6.95 (T) 

-6.33 (E) 

0.62 

0.89 (X) 

-7.63 (P) 

-5.69 (E) 

1.94 

2.38 (E) 

2H-TiSe2 

-6.31 (P) 

-5.89 (E) 

0.42 

0.77 (E) 

-6.63 (P) 

-5.15 (M) 

1.48 

2.13 (M) 

2H-TiTe2 

-5.40 (P) 

-5.44 (E) 

0.00 

0.31 (E) 

-5.52 (T) 

-5.06 (E) 

0.45 

1.21 (T) 

2 H-WO 2 

-6.73 (P) 

-5.41 (K) 

1.32 

1.65 (P) 

-7.38 (P) 

-4.73 (K) 

2.65 

3.18 (P) 

2 H-WS 2 

-5.75 (K) 

-4.24 (K) 

1.51 

1.51 (K) 

-6.28 (K) 

-3.85 (K) 

2.43 

2.43 (K) 

2H-WSe2 

-5.13 (K) 

-3.91 (K) 

1.22 

1.22 (K) 

-5.61 (K) 

-3.53 (K) 

2.08 

2.08 (K) 

2H-Zr02 

-8.44 (T) 

-6.85 (T) 

1.59 

1.70 (T) 

-9.71 (T) 

-5.63 (T) 

4.08 

4.19 (T) 

1T-Zr02 

-8.20 (T) 

-3.82 (K) 

4.37 

4.63 (P) 

-9.73 (X) 

-1.97 (M) 

7.76 

8.25 (E) 

2H-ZrS2 

-7.02 (T) 

-6.18 (X) 

0.85 

1.03 (X) 

-8.02 (T) 

-5.56 (E) 

2.46 

2.69 (T) 

lT-ZrS 2 

-6.58 (P) 

-5.55 (M) 

1.03 

1.53 (P) 

-7.60 (P) 

-4.72 (E) 

2.88 

3.61 (P) 

2H-ZrSe2 

-6.47 (T) 

-5.82 (X) 

0.64 

0.91 (X) 

-7.29 (T) 

-5.33 (E) 

1.96 

2.27 (X) 

lT-ZrSe 2 

-5.66 (P) 

-5.41 (M) 

0.25 

0.87 (P) 

-6.53 (P) 

-4.68 (M) 

1.85 

2.63 (P) 

2H-ZrTe2 

-5.62 (P) 

-5.44 (E) 

0.18 

0.47 (E) 

-6.17 (T) 

-5.16(E) 

1.01 

1.41 (E) 


“See Fig. lb 
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Table 4: LDA band gaps (in eV) as function of strain. The character in the paranthesis denotes 
whether the gap is (I)ndirect or (D)irect. 


name 

-2% 

-1% 

0 % 

1 % 

2 % 

2 H-Cr02 

0.77 (I) 

0.59 (1) 

0.42 (I) 

0.28 (I) 

0.15(1) 

2 H-CrS2 

1.04(D) 

0.98 (D) 

0.92 (D) 

0.86 (I) 

0.72 (I) 

2 H-CrSe2 

0.84 (D) 

0.79 (D) 

0.74 (D) 

0.70 (D) 

0.67 (D) 

2 H-CrTe2 

0.59 (D) 

0.56 (D) 

0.52 (D) 

0.49 (D) 

0.47 (D) 

2 H-Ge02 

1.77 (I) 

1.57 (1) 

1.37 (I) 

1.19(1) 

1.00 (I) 

2H-Hf02 

2.02 (I) 

1.96 (1) 

1.89 (I) 

1.82 (I) 

1.75 (I) 

2H-HfS2 

0.96 (I) 

0.96 (1) 

0.96 (I) 

0.94 (I) 

0.93 (I) 

2 H-HfSe2 

0.71 (I) 

0.79 (1) 

0.80 (I) 

0.80 (I) 

0.79 (I) 

2 H-Hfre2 

0 .10(1) 

0.21 (1) 

0.31 (I) 

0.41 (I) 

0.50 (I) 

2 H-M 0 O 2 

1.41 (I) 

1.15(1) 

0.91 (I) 

0.70 (I) 

0.50 (I) 

2 H-M 0 S 2 

1.82 (I) 

1.78 (D) 

1.65 (I) 

1.41 (I) 

1.19(1) 

2 H-MoSe2 

1.50 (I) 

1.53 (D) 

1.44 (D) 

1.34 (D) 

1.26(D) 

2 H-MoTe2 

1.16(1) 

1.16(D) 

1.07 (D) 

1.00 (D) 

0.93 (D) 

2H-Sc02 

1 .12(1) 

1.13(1) 

1.15(1) 

1.16(1) 

1.17(1) 

2H-ScS2 

0.48 (I) 

0.50 (1) 

0.50 (I) 

0.50 (I) 

0.49 (I) 

2 H-ScSe2 

0.31 (I) 

0.34 (1) 

0.36 (I) 

0.38 (I) 

0.37 (I) 

2 H-Sn02 

0.83 (I) 

0.69 (1) 

0.57 (I) 

0.43 (I) 

0.31 (I) 

2 H-SnS2 

0.56 (I) 

0.60 (1) 

0.63 (I) 

0.64 (I) 

0.65 (I) 

2H-Ti02 

1.25 (I) 

1.19(1) 

1.14(1) 

1.07 (I) 

1.00 (I) 

2H-TiS2 

0.64 (I) 

0.64 (1) 

0.63 (I) 

0.61 (I) 

0.58 (I) 

2 H-TiSe2 

0.33 (I) 

0.42 (1) 

0.51 (I) 

0.53 (I) 

0.52 (I) 

2 H-TiTe2 

M 

M 

0.09 (I) 

0.17 (I) 

0.24 (I) 

2 H-VSe2 

M 

M 

M 

M 

M 

2 H-VTe2 

M 

M 

M 

0.11 (I) 

0.17 (I) 

2 H-WO 2 

1.93 (I) 

1.63 (1) 

1.36(1) 

1 .10(1) 

0.87 (I) 

2 H-WS 2 

1.88 (I) 

1.94(1) 

1.80 (D) 

1.58 (I) 

1.34 (I) 

2 H-WSe2 

1.55 (I) 

1.61(1) 

1.54 (D) 

1.43 (D) 

1.32(D) 

2 H-Zr02 

1.76 (I) 

1.70 (1) 

1.63 (I) 

1.55 (I) 

1.47 (I) 

2 H-ZrS2 

0.87 (I) 

0.87 (1) 

0.86 (I) 

0.84 (I) 

0.82 (I) 

2 H-ZrSe2 

0.69 (I) 

0.70 (1) 

0.70 (I) 

0.70 (I) 

0.69 (I) 

2 H-ZrTe2 

0.18 (I) 

0.29 (1) 

0.38 (I) 

0.44 (I) 

0.45 (I) 

1 T-Ge02 

4.07 (I) 

3.81 (1) 

3.56 (I) 

3.31 (I) 

3.06 (I) 

lT-GeS 2 

0.48 (I) 

0.53 (1) 

0.57 (I) 

0.61 (I) 

0.64 (I) 

1T-Hf02 

4.63 (I) 

4.65 (1) 

4.66 (I) 

4.65 (I) 

4.56 (I) 

lT-HfS2 

0.90 (I) 

1.04(1) 

1.16(1) 

1.27 (I) 

1.38 (I) 

lT-HfSe 2 

0.26 (I) 

0.40 (1) 

0.52 (I) 

0.64 (I) 

0.75 (I) 

1 T-Mn02 

0.69 (I) 

0.70 (1) 

0.72 (I) 

0.75 (I) 

0.77 (I) 

1T-Ni02 

1.30 (I) 

1.24(1) 

1.18(1) 

1.11 (I) 

1.05 (I) 

lT-NiS2 

0.35 (I) 

0.48 (1) 

0.55 (I) 

0.60 (I) 

0.65 (I) 

lT-NiSe 2 

M 

M 

0.16(1) 

0.23 (I) 

0.29 (I) 

1T-Pb02 

1.55 (I) 

1.44(1) 

1.32(1) 

1.20 (I) 

1.08 (I) 

lT-PbS2 

0.56 (I) 

0.61 (1) 

0.64 (I) 

0.67 (I) 

0.69 (I) 

1T-Pd02 

1.46 (I) 

1.39 (1) 

1.32(1) 

1.23 (I) 

1.15(1) 

lT-PdS2 

1.09 (I) 

1.13(1) 

1.17(1) 

1.14(1) 

1.06 (I) 

lT-PdSe 2 

0.55 (I) 

0.61 (1) 

0.66 (I) 

0.71 (I) 

0.72 (I) 

lT-PdTe 2 

M 

0.15 (1) 

0.22 (I) 

0.27 (I) 

0.32 (I) 

1T-Pt02 

1.78 (I) 

1.68 (1) 

1.59(1) 

1.50 (I) 

1.41 (I) 

lT-PtS2 

1.73 (I) 

1.71(1) 

1.66 (I) 

1.61 (I) 

1.54 (I) 

lT-PtSe 2 

1.20 (I) 

1.25 (1) 

1.29 (I) 

1.25 (I) 

1.17(1) 

lT-PtTe 2 

0.50 (I) 

0.63 (1) 

0.69 (I) 

0.74 (I) 

0.73 (I) 

1 T-Sn02 

2.89 (I) 

2.72 (1) 

2.54 (I) 

2.36 (I) 

2.18 (I) 

lT-SnS 2 

1.35 (I) 

1.38(1) 

1.41 (I) 

1.43 (I) 

1.45 (I) 

lT-SnSe2 

0.58 (I) 

0.63 (1) 

0.67 (I) 

0.71 (I) 

0.74 (I) 

1T-Ti02 

2.82 (I) 

2.74 (1) 

2.66 (I) 

2.58 (I) 

2.50 (I) 

1 T-Zr02 

4.48 (I) 

4.50 (1) 

4.38 (I) 

4.23 (I) 

4.10(1) 

lT-ZrS 2 

0.84 (I) 

0.96 (1) 

1.08 (I) 

1.19(1) 

1.29 (I) 

lT-ZrSe 2 

0.17 (I) 

0.30 (1) 

0.42 (I) 

0.53 (I) 

0.64 (I) 
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Table 5: Effective electron and hole masses together with the direct and indirect exciton masses 
calculated from the GqWo quasiparticle band structures with spin-orbit interaction included". The 
slope of the quasi-2D dielectric function at <? = 0 is is shown and the exciton binding energies ob¬ 
tained from a quasi-2D Mott-Wannier model. 


name 

m* (me) 

m* (me) 

Met“‘ (me) 

^indirect 

de^ldq\^^Q 

(A) 

^direct (gY) 

^.indirect 

2H-Cr02 

1.1/1.1 

1.3/1.3 

0.6/0.6 

0.59/0.59 

45 

0.56 

0.56 

2H-CrS2 

1.1/1.1 

0.88/0.88 

0.49/0.49 

0.49/0.49 

61.5 

0.42 

0.42 

2H-CrSe2 

1.1/1.1 

0.97/0.97 

0.52/0.52 

0.52/0.52 

73.6 

0.37 

0.37 

2H-CrTe2 

0.99/0.99 

0.89/0.89 

0.47/0.47 

0.47/0.47 

95.4 

0.30 

0.30 

2H-Ge02 

0.32/0.32 

5/8.6 

0.3/0.3 

0.3/0.31 

10.2 

1.24 

1.25 

1T-Ge02 

0.34/0.34 

4.1/2.6 

0.81/0.81 

0.32/0.3 

7.39 

1.97 

1.55 

lT-GeS 2 

0.64/0.22 

1.5/0.25 

0.82/0.13 

0.26/0.14 

27.6 

- 

- 

2H-Hf02 

2.5/1.7 

4.7/4.5 

1.1/1.5 

1.6/1.2 

9.85 

- 

- 

1T-Hf02 

3.3/0.79 

1.1/3.1 

1.1/0.79 

0.53/1.1 

8.82 

- 

- 

2H-HfS2 

-21/1.2 

3/0.95 

1.5/4.5 

2.4/0.56 

23.5 

- 

- 

lT-HfS2 

1.4/0.29 

0.63/0.63 

0.78/0.78 

0.44/0.2 

27.7 

0.80 

- 

2H-HfSe2 

-38/0.75 

2.9/0.51 

0.97/4.7 

1.9/0.33 

32.6 

- 

- 

lT-HfSe 2 

1.8/0.23 

0.51/0.51 

0.47/0.47 

0.4/0.16 

41.6 

0.55 

- 

2H-Hfre2 

-48/0.46 

0.93/1.4 

-4.2/0.23 

1.1/0.33 

66.7 

- 

- 

2 H-M 0 O 2 

0.51/0.51 

0.8/0.8 

0.75/0.75 

0.31/0.31 

31.4 

0.75 

0.62 

2 H-M 0 S 2 

0.55/0.55 

0.56/0.56 

0.28/0.28 

0.28/0.28 

44.3 

0.47 

0.47 

2H-MoSe2 

0.49/0.49 

0.61/0.61 

0.27/0.27 

0.27/0.27 

51.2 

0.42 

0.42 

2H-MoTe2 

0.65/1.1 

0.64/0.64 

0.31/0.31 

0.32/0.4 

65.4 

0.36 

- 

1T-Ni02 

1.1/2.1 

4.2/33 

0.62/1.4 

0.87/1.8 

35.8 

- 

- 

lT-NiS2 

0.39/0.79 

1.3/1.4 

1.4/0.35 

0.3/0.51 

79.3 

- 

- 

lT-NiSe 2 

0.29/0.52 

-58/2.9 

0.26/0.98 

0.29/0.44 

121 

- 

- 

1T-Pb02 

0.39/0.39 

53/5 

0.41/0.41 

0.39/0.36 

12.8 

1.20 

1.17 

lT-PbS2 

0.83/0.33 

10/0.64 

0.96/0.2 

0.48/0.27 

32 

- 

- 

1T-Pd02 

1.3/2.6 

6.8/1.1 X 10^ 

0.58/1.4 

1.1/2.5 

28.1 

- 

- 

lT-PdS2 

0.35/0.77 

1.3/0.74 

0.29/1.7 

0.26/0.41 

54.4 

- 

- 

lT-PdSe 2 

0.28/0.52 

6.5/7.1 

0.22/1.3 

0.27/0.49 

76.1 

- 

- 

lT-PdTe 2 

0.23/0.31 

0.99/0.99 

0.16/0.31 

0.19/0.24 

134 

- 

- 

1T-Pt02 

1.1/2.2 

1.6/16 

0.47/1.1 

0.65/2 

21.2 

- 

- 

lT-PtS2 

0.32/0.8 

0.48/2.1 

0.27/0.62 

0.19/0.58 

38.7 

- 

- 

lT-PtSe 2 

0.26/0.56 

1.1/0.6 

0.25/6.6 

0.2/0.31 

50.2 

- 

- 

lT-PtTe 2 

0.23/0.38 

1.6/1.6 

0.19/0.42 

0.2/0.31 

75.3 

- 

- 

2H-Sn02 

0.31/0.31 

6.3/11 

0.33/0.33 

0.29/0.3 

12 

1.16 

1.13 

1T-Sn02 

0.33/0.33 

3/4.5 

0.55/0.82 

0.3/0.31 

8.11 

- 

1.45 

2H-SnS2 

0.69/0.34 

2.3/2.3 

0.2/0.2 

0.53/0.29 

24.2 

0.64 

- 

lT-SnS 2 

0.74/0.28 

2.8/0.32 

0.9/0.16 

0.34/0.18 

21.6 

- 

- 

lT-SnSe 2 

0.67/0.24 

2.2/0.26 

0.87/0.14 

0.29/0.15 

31.8 

- 

- 

2H-Ti02 

1.4/2.2 

5/3.6 

1.5/1.9 

1.1/1.4 

14.2 

- 

- 

1T-Ti02 

8.2/4.6 

1.1/4.8 

1.3/1.1 

0.99/2.4 

14.5 

- 

- 

2H-TiS2 

-29/0.78 

1/1 

-21/0.83 

1.1/0.45 

38.1 

- 

- 

2H-TiSe2 

8.1/0.52 

0.63/0.63 

47/0.32 

0.59/0.29 

55.6 

- 

- 

2H-TiTe2 

3.4/0.4 

1.1/0.69 

1.6/0.83 

0.74/0.26 

129 

- 

- 

2 H-WO 2 

0.45/0.45 

0.76/0.76 

0.78/0.78 

0.28/0.28 

26.7 

0.84 

0.68 

2 H-WS 2 

0.46/0.46 

0.42/0.42 

0.22/0.22 

0.22/0.22 

39.9 

0.48 

0.48 

2H-WSe2 

0.48/0.48 

0.44/0.44 

0.23/0.23 

0.23/0.23 

46.2 

0.43 

0.43 

2H-Zr02 

2.7/1.1 

4.1/4.2 

1 .1/1 

1.6/0.87 

11 

1.59 

- 

1T-Zr02 

3.6/I.4 

1.1/3.7 

1.2/1.1 

0.73/1.3 

9.63 

1.82 

- 

2H-ZrS2 

1.4/3.1 

2.6/0.95 

0.58/3 

1.1/0.65 

24.9 

- 

- 

lT-ZrS 2 

2/0.34 

0.71/0.71 

0.67/0.67 

0.53/0.23 

30.6 

0.73 

- 

2H-ZrSe2 

1.3/1.4 

2.2/0.59 

0.51/2 

0.82/0.41 

34.1 

- 

- 

lT-ZrSe 2 

1.9/0.28 

0.59/0.59 

0.47/0.47 

0.45/0.19 

48 

0.50 

- 

2H-ZrTe2 

2.1/0.53 

1.8/2 

3.4/0.13 

0.99/0.41 

74.9 

- 

- 


“Negative masses occur in some directions due to bad fitting. This is usually the case if the band structure is very flat in one 


direction but highly varying in the other direction. Thus negative masses generally mean that the mass in this direction is 

much larger than in the other direction. 
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